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o 
^. 

Mh, Using supersymmetric quantum mechanics we construct the quasi- 

-^ I exactly solvable (QES) potentials with arbitrary two known eigen- 

j^ ■ states. The QES potential and the wave functions of the two energy 

^ . levels are expressed by some generating function the properties of 

^ ' which determine the state numbers of these levels. Choosing different 

generating functions we present a few explicit examples of the QES 
S^ . potentials. 
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1 Introduction 

From the earlier days of quantum mechanics there has been continual inter- 
est in the models for which the corresponding Schrodinger equation can be 
solved exactly. The number of totaly exactly solvable potentials is rather lim- 
ited. Therefore, recently much attention has been given to the quasi-exactly 
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solvable (QES) potentials for which a finite number of the energy levels and 
the corresponding wave functions are known in the explicit form. 

The first examples of QES potentials were given in [^ ^ 0, ^. Subse- 
quently several methods were worked out for generating QES potentials and 
as a result many QES potentials were established p|, ||, |^, ||, ||, |l^, [Tl|, |l^, |l^]. 
One of the methods is the generation of new QES potentials using supersym- 
metric (SUSY) quantum mechanics [l^, 0, ^ (for review of SUSY quantum 
mechanics see |jl6|). The idea of the SUSY method for constructing QES po- 
tentials is the following. Starting from some initial QES potential with n+1 
known eigenstates and using the properties of the unbroken SUSY one ob- 
tains the supersymmetric partner potential which is a new QES potential 
with the n known eigenstates. 



In our recent paper [|T^ we have proposed a new SUSY method for gener- 
ating QES potentials with two known eigenstates. This method, in contrast 
to the one in papers |]13|, |I^, |T^, does not require the knowledge of the initial 
QES potential for generating a new QES one. The general expression for 
the superpotential, the potential energy and two wave functions which corre- 
spond to two energy levels were obtained. Within the frames of this method 
we have obtained QES potentials for which we have found in the explicit form 
the energy levels and wave functions of the ground and first excited states. 
One should mention here also paper [|I8| where the general expression for the 
QES potentials with two known eigenstates was obtained without resorting 
to the SUSY quantum mechanics (see also a very recent paper by Dolya and 
Zaslavskii |19]). Although this method is direct and simpler than the SUSY 



approach the latter still has some advantages. Namely, the SUSY method 



developed in |T^ can be extended for the generation of QES potentials with 
three known eigenstates [^ and conditionally exactly solvable potentials 
(CES) [^. The CES potentials are those for which the eigenvalue problem 
for the corresponding Hamiltonian is exactly solvable only when the poten- 
tial parameters obey certain conditions [^. About using SUSY quantum 
mechanics for the construction of the CES potentials see 
worth to mention the very recent paper 



23, 24, m. It is 



where the authors established 
the connection between SUSY approach for constructing QES potentials with 
two known eigenstates 



17| , ^^ and the Turbiner approach 



Note that the general expressions for QES potential and corresponding 
wave functions derived by Dolya and Zaslavskii in [|T9l without resorting 



to SUSY quantum mechanics are the same as those obtained in O, ET 



using SUSY method. In |T^, ^ we have used these general expression for 
constructing QES potentials with the ground and first excited states. A new 
interesting result obtained by Dolya and Zaslavskii is that they have shown 
that it is possible to obtain not only the ground and first excited states but 
any pair of the energy levels and the corresponding wave functions. 

The aim of the present paper is to extend the SUSY method proposed in 



our papers |17, 21| for constructing QES potentials with arbitrary two known 



eigenstates. In JT^ ^ we used nonsingular superpotentials and obtained 



QES potentials with explicitly known ground and first excited states. In 
the present paper using singular superpotentials we obtain nonsingular QES 
potentials for which we know in the explicit form any pair of the energy 
levels and the corresponding wave functions. Of course, the idea of using 



singular superpotential in SUSY quantum mechanics is well know (see [|T6 
for review). Nevertheless, a new moment of present paper is that in the 
frame of SUSY quantum mechanics with singular superpotentials we derive 
nonsingular QES potentials with arbitrary two known eigenstates. 

2 SUSY quantum mechanics and QES prob- 
lems 

In the Witten's model of supersymmetric quantum mechanics the SUSY 
partner Hamiltonians H± read 



where 



H± = B^B^ = -\^,+V^{x\ 


(1) 


\ 1 d \ 
B^= '{t, +Wix)], 


(2) 


l{w\.)±w'i.)), w'i.)-''^/\ 


(3) 



V±{x) 

and W{x) is referred to as a superpotential. In this paper we shall consider 
the systems on the full real line — cxo < x < oo. 

We shall study the eigenvalue problem for the Hamiltonian if_ 

B+B-ij^{x)=E^P^{x). (4) 



Applying to the left and right hand sides of this equation the operator B 
we obtain the equation for the eigenvalue problem of the Hamiltonian H^ 

B-B+{B-rE{x)) = E{B'rE{x)) (5) 

from which follows that 

V'+(x) = CB-rE{x) (6) 

is the solution of the eigenvalue problem for the Hamiltonian if+ with the 
energy E, C = l/^/E is the normalisation constant if the wave function 
'il)~^{x) is square integrable. It is also possible that 4'^{x) does not satisfy 
the necessary conditions and does not belong to the eigenf unctions of the 
Hamiltonian H^. Such situation takes place for singular superpotentials (see 



review [|l6l). Nevertheless, the function i{j~^{x) is the solution of the equation 



(H). Applying to equation (||) the operator 5+ we obtain 

^^(a;)=Cfi+^+(a;). (7) 

Now let us analyse the eigenvalue problem for the Hamiltonian H_. Due 
to the factorization of the Hamiltonian the wave function of the zero energy 
state satisfies the equation B~ipQ = and reads 

^o-(x) = Co exp (^- J W{x)dx^ , (8) 

Cq is the normalization constant. In order to satisfy the condition of the 
square integrability of the wave function (^ we put 

sign(Vr(±cx))) = ±1. (9) 

We are interested in the potential energy VI (x) free of singularities. The 
simplest way to satisfy this condition is to consider a superpotential W{x) 
which is free of singularities. Then iPq{x) corresponds to the zero energy 
ground state of the Hamiltonian H_. Just this case was considered in our 
paper fl^. But it is also possible to get a nonsingular potential energy 
V-{x) using a singular superpotential. Let us assume that W{x) has the 
simple poles at the points Xk with the following behaviour in the vicinity of 

Xk 

W{x) = ^-^ + Ao + Ai{x- Xk) + 0{{x - Xkf). (10) 



Then V- (x) in the vicinity of Xk reads 
2V.{x) = ^-/^;^-^ + ^^ + 2^^ + 2A.,A, + Al-A,+ 0{x - x,). (11) 

yX Xi^j X Xj^ 

The first case A_i = leads to the nonsingular W{x) and V-{x). The second 
case A_i = —1 and Aq = gives the nonsingular potential energy V-{x) with 
a singular superpotential. Here it is worth stressing that the SUSY partner 
potential energy in this case is singular with the following behaviour in the 
vicinity of Xk'- 2V+{x) = 2/(x — x^)^ — Ai + 0{x — Xk). 

Let us analyse the second case. Using equation (H) we obtain the be- 
haviour of the wave function in the vicinity of the point x^ as follows ipQ (x) ~ 
\x — Xk\{l — Ai{x — XkY/2). As we see d'ip^{x)/dx is a discontinuous func- 
tion at the points Xk- In order to obtain the wave function with continuous 
derivative we use the simple fact that if '?/'o"(a^) in the domain Xk < x < Xk+i 
satisfies the Schrodinger equation then the function with the opposite sign 
—xl^Q^x) satisfies the same equation too. Thus, we can change the sign of 
wave function in the domains {xk,Xk+i) in such a way that iPq{x) and its 
derivative dx(:Q{x)/dx will be continuous functions. In fact it means that it 
is necessary to make the substitution |/| -^ f. Then the behaviour of the 
wave function in the vicinity of x^ is 

^o~(x) ~ (x - Xk){l - A,{x - Xkf/2) (12) 

and the wave function has zeros at the points Xk- Thus, the zero energy 
wave function has n nodes (n is the number of poles of the superpotential) 
and corresponds to the n-th excited state. Note that in this case the ground 
state energy is less than zero. Choosing different superpotentials W{x) we 
can easily construct different QES potentials with one known eigenstate. 

In contrast to this the construction of the QES potentials with two known 
eigenstates is not a trivial problem. One state of the Hamiltonian H„ with 
the zero energy is known and is given by (§). In order to obtain one more 
state of iJ_ we use the following well-known procedure exploited in SUSY 
quantum mechanics. Let us consider the SUSY partner of if_, i.e. the 
Hamiltonian if+. If we calculate some state of if+ we immediately find a 
new excited state of i7_ using transformation (0). In order to calculate some 
state of if+ let us rewrite it in the following form 

H+ = H^y + e = B+B^ + e, e > 0, (13) 



which leads to the following relation between the potential energies 

V+{x) = Vl^\x) + e, (14) 

and superpotentials 

W\x) + W'{x) = W^{x) - W[{x) + 2e, (15) 

where e is the energy of the state of iJ+ since we supposed that H_ similarly 
to H_ has the zero energy state, Bf and V_ (x) are given by (@) and (|^) 
with the new superpotential Wi{x). 

As we see from (|13D the wave function of H^ with the energy E = e 
is also the zero-energy wave function of H_ and it satisfies the equation 
Biip'^{x) = 0. The solution of this equation is 

^+(x) = C+ exp (- f Wi{x)dx) . (16) 

Using (0) we obtain the wave function of the excited state with the energy 
level E = e for the Hamiltonian H_ 

^-{x) = C- W+{x) exp (- f Wi{x)dx) (17) 

where we have introduced the notation W+{x) = Wi{x) + W{x). In order 
to satisfy the square integrability of this function at infinity we impose the 
superpotential Wi{x) with the same condition as W{x) @. Then W+{x) 
satisfies the same condition (P) too. 

In order to obtain the explicit expression for the wave function iPq{x) 
with zero energy and the wave function ip' (x) with the energy e given by (|^) 
and ( p!?!) it is necessary to obtain the explicit expression for superpotentials 
W{x) and Wi{x). This is the subject of the next section. 

3 Solutions for superpotentials and construc- 
tion of nonsingular QES potentials 

The superpotentials W{x) and Wi{x) satisfy equation ([151) . Note, that (|T5|) 
is the Riccati equation which cannot be solved exactly with respect to Wi{x) 
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for a given W{x) and vice versa. But we can find such a pair of W{x) and 



Wi{x) that satisfies equation (plsl). For this purpose let us rewrite equation 
(|T3p in the following form 

iy|(x) = W4x)W+{x) + 2e, (18) 

where 

W+{x) = Wi{x) + W{x), (19) 

W4x)=Wi{x)-W{x). (20) 

This new equation can be easily solved with respect to W- (x) for a given 
W + (x) and vice versa. In this paper we use the solution of equation (|l^) 
with respect to W^{x) 

W^x) = (W^ix) - 2e)/W+{x). (21) 

Then from {^^, (|l^) and (^) we obtain the pair of W{x), Wi{x) that satisfies 
equation (|T^) 

W{x) = i {w+{x) - «(a;) - 2e)/W+ix)) , (22) 

Wiix) = ^ (w+ix) + iW[ix) - 2e)/W+ix)) , (23) 

here W+(x) is some function of x which generates the superpotentials W{x) 
and Wi{x). Let us stress that W{x), Wi{x) and W+{x) must satisfy con- 
dition (H). It is necessary to note that the general solutions (0), (p3|) of 



equation (^) was obtained earlier in [^] in the context of parasupersym- 
metric quantum mechanics. 

In our earlier paper [jl^ we considered only free of the singularities su- 
perpotentials W{x) and Wi{x). To satisfy a nonsingularity of the superpo- 
tentials we considered a continuous function W+{x) that has only one simple 
zero. Because W+{x) was considered as continuous function which satisfies 
condition (|]) the function W+{x) must have at least one zero. Then, as we 
see from (|2l|), (||) and (|3D, W_{x), W{x) and Wi{x) have the poles. In 



order to construct the superpotentials free of singularities we supposed that 
W^{x) has only one simple zero at a: = Xq. In this case the pole of W^{x) 
and W{x), Wi{x) at x = xq can be cancelled by choosing e = PV|(xo)/2. 
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Now let us consider more general cases of the function W+{x) which leads 
to nonsingular QES potential energy V-{x). 

Case 1. Suppose that W+{x) has simple zeros at the points x^, /c = 1, ..., n 

W+{x) = W'_^{xk){x - Xk) + \w'l{xk){x - Xkf + 0{{x - Xkf). (24) 

The zeros of the function W-^{x) lead to the poles of the function W-{x) and 
the superpotential W{x) with the following behaviour in the vicinity of Xk 

(25) 
The behaviour of the superpotential Wi{x) in the vicinity of Xk is similar to 
W{x) only with the opposite sign. It is worth comparing the superpotential 
(II) with (|TOD: A_i = e/W'+{xk) - 1/2 and A^ = -W"{xk){e/W'^{xk) + 
l/2)/2W'{xk)- This superpotential leads to the following behaviour of the 
potential energy in the vicinity of x^ 



2VJx] 



1 



W'+{xk)) 4 



^"^^^) ' ^■O(const). 



(X - XkY W'{Xk) {X - Xk) , 

(26) 
Thus in the case 

W'^{xk) = ±2e (27) 

the potential energy V-{x) is free of singularities. It is convenient to divide 
the set of Xk into two subsets x^ {k = 1, ...,n+) and x^ {k = 1, ...,n~) for 
which W^(xj^) = 2e > and iy|(x^) = — 2e < 0. We suppose in this paper 
that e > 0. Because of W|(x^) = 2e the singularity at the points x'l is 
cancelled and W{x), Wi{x) have singularities only at the points x^ 

W{x)= ~^_ +0(x-Xfc), (28) 

Wi{x) = — ^ + 0(x-Xfc). (29) 

X — Xj. 

Substituting W(x) into (|) and using the result of the previous section (see 
Eq. ([T^) ) we see that the wave function V'o^l^) with zero energy has n~ zeros 
at the points x^, namely i/jq ~ (x — x^) in the vicinity of x^. Substituting 



W+{x) given by (^) and Wi{x) given by (^) into (|17|) we obtain that the 
wave function ip' (x) with the energy e has n+ zeros at the points x^: ip' (x) ~ 
(x — x'l). When W+{x) is the continuous function satisfying condition (|^) 
then n^ = n~ + 1. Thus in this case ipoix) and '?/'7(a;) correspond to n~-th 
and (n~ + l)-th excited states respectively. 

Case 2. Now let us assume that the function W+{x) in addition to the 
zeros has the simple poles at the points a;° with the behaviour in the vicinity 
of a;° as follows 

W+{x) = -^ + Go + 0{x - xl). (30) 






Then 



W{x) = -.—^ + :^77^(G-i - 1) + 0{x - xl), (31) 

1 G^i -^ , 1 Go 

IcT 



W^i(^) = o— ^ + 97^(^-1 + 1) + ^(^ - ^°)' (32) 

.Z X Xt, .Z Or 1 



here we drop out the terms of the order [x — a;°). Note that Go and G_i can 
depend on k. For the sake of simplicity we are omitting this dependence. 
Comparing superpotential (|3T|) and (|^) we conclude that the superpotential 
W{x) ( pT] ) gives a nonsingular potential energy for the case 2a: G_i = — 1 
and Go is an arbitrary constant, and for the case 2b: G_i = —3 and Gq = 0. 
Case 2a. For the case G_i = —1 in the vicinity of x° we have nonsingular 
superpotential W{x) and singular Wi{x) 

W{x)=Gq + 0{x-xI), (33) 

Wi{x) = ^^ + 0{x-xl). (34) 

X X j^ 

The wave functions iIjq^x) and ip~{x) calculated with these superpotentials 
do not have the zeros at the points x°, which we will denote in this case as 
Ofc, fc = 1, ...,n°. Nevertheless, the important moment is that now Wj^{x) in 
additional to n = ri^ + n~ zeros at points x~l and x^ has n^ poles at the 
points x° and thus is not a continuous function. As a result in this case we 
have n"*" = n~ + n° + 1. Using the result obtained in the case 1 we see that 
iIjq{x) has the zeros at x^ and corresponds to the n~-th excited state and 
iIj~ has the zeros at x~l and corresponds to the {n~ +n° + l)-th excited state. 



Case 2b. The case G_i = —3 and Go = leads to the following behaviour 
of superpotentials in the vicinity of x° 

Wix) = ^^ + Oix-xl), (35) 

k 

Wi{x) = ^^ + 0(x - xl) (36) 

The wave functions ipQ^x) and ip'^x) calculated with these superpotentials 
have common zeros at the points x^, which we will denote in this case as 
bk, k = l,...,mP. Thus, when in addition to poles the function W+{x) has 
n = n^ + n~ zeros at the the points x^ and x^ the wave function iPq{x) 
corresponds to {n~ + m^)-th excited state and ip'ix) corresponds to {n~ + 
2m° + l)-th excited state. 

Let us consider the general case which combines the cases 1, 2a and 2b. 
The function W^{x) has n~ zeros with negative derivatives at the points 
x^ {k = 1, ...,n'^), n^ poles at the points a^ {k = 1, ...,n^) with asymptotic 
behaviour in the vicinity of these points —l/{x — a^) + const and m^ poles 
at the points bk {k = 1, ...,m°) with the asymptotic behaviour —S/{x — bk) ■ 
The number of zeros of the function Wj^{x) with a positive derivative at the 
points x~l is the following n+ = n~ + n° + mP + 1. The wave function iPq{x) 
has the (n~ + mP) nodes at the points x~^ and bk and thus it corresponds to 
the {n^ + m'')-th excited state. The wave function ip~{x) has n^ + ivP = 
n~ +n^ + 2m^ + 1 nodes at the points x^ and bk and thus it corresponds to 
the (n~ + n° + 2m^ + l)-th excited state. 

Thus, the considered cases 1, 2a, 2b and the combined general case lead 
to the nonsingular QES potential energy V^{x) given by (^, where the su- 
perpotential W{x) is expressed over the function iy+(x) by equation (|2^) . 
The zero energy wave function ip'^{x) and the wave function 'ip^{x) with the 
energy e are given by (H) and (0), respectively. Note, that in the case of non- 
singular superpotential W{x) the zero energy wave function corresponds to 
the ground state, but in the case of singular superpotential the energy of the 
ground state is less than zero and the zero energy wave function corresponds 
to an excited state. 

In the conclusion of this section let us discuss the second possibility for 
construction of the QES potentials with two known eigenstates. Namely, we 
can use the solution of equation (|18D with respect to W^{x) [^. We found 
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that 

^ ' 0'(x)' ^^ ' (f)'{x) ^ ' 

satisfies equation ([18|) , where 0(x) is new generating function. Using VF+ and 
W^_ given by (^) we obtained the QES potential and two wave functions in 
term of 0(x). Choosing generating functions 0(x) with one zero we obtained 
QES potentials with explicitly know ground and first excited states. Note 
that basic equation derived by Dolya and Zaslavskii in [^ without resorting 



to SUSY quantum mechanics are the same as was earlier obtained in pT 



using SUSY method (in [|T9| 0(x) is denoted as ^(x)). A new result obtained 
by Dolya and Zaslavskii is that they have shown how one can obtain not only 
the ground and first excited states but any pair of states using generating 
function 0(a;) and their derivative (j)'{x) with zeros and poles. In the present 
paper we are working with generating function Wj^{x) which is related with 
cf)[x) by (^). As we see zeros and poles of 0(x) (or ^(x)) lead the zeros of 
Wj^{x) and zeros of 0'(x) (or i'{x)) lead to the poles of W^{x). 

4 Examples of QES potentials 

Note, that all expressions depend on the function Wj^{x). We may choose 
various functions W+{x) and obtain as a result various QES potentials. Note 
also that when the function Wj^{x) generates the potential energy V-{x) 
then W+{x/a)/a generates the potential energy V_{x/a)/a?. This scaling is 
useful for comparing, in principle, the same potential energies but written in 
different forms as a result of a different measurement units. 

To illustrate the above described method we give two explicit examples 
of the nonsingular QES potentials. 

Example 1 

Let us consider a continuous function Wj^{x) which corresponds to case 

1 

x^ - 1 
W+{x) = oix^—- (38) 

x"' + 1 

where a > 0. This function has three zeros at the points 0,±1. The de- 
nominator is written in order to satisfy condition (p7|), namely W^(0) = —a, 
M/^|(±l) = a. Note that n+ = 2, n~ = 1 and thus we have QES potential 
with explicitly known first and second excited states. From the condition of 
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nonsingularity of the potential energy it follows that e = a/2. Then using 
( p2[) and (|23| ) we obtain for the superpotentials 

W{x) = '^x + {l-a)-^--, (39) 

1 X^ + 1 X 

ry X \ 

iyi(a:) = -a:-(l + a)^^^ + -. (40) 

The superpotential W{x) gives the following QES potential 

2V_(x) = ^x2-(l-a)(3-a)^-^l^-2a(l-a)-^ + a(l-a)-^a. 

(41) 
The zero energy wave function (^ and wave function with energy e = a/2 
(|13) read 

^-(x) = Cox(x2 + l)("-i)/2 exp (-^x^) , (42) 

^-[x) = C,{x^ - l){x^ + l)("-i)/2g^p ("-^x^") . (43) 

As we see iPq{x) has one node and thus really corresponds to the first excited 
state, ip~{x) has two nodes and corresponds to the second excited state. 

Note that QES potential (^) has similar structure as the potential stud- 
ied in [|19| but in fact it is another potential. In our case we have the QES 
potential with explicitly known first and second eigenstates, whereas for the 
QES potential studied in [1^ the ground and second excited states are explic- 



itly known. It is interesting to note also that in our case QES potential (p]) 
at the value of parameter a = 1 becomes exactly solvable and corresponds 
to the harmonic oscillator. 

Example 2 

Let us consider a more complicated example for which the function Wj^{x) 
has two poles and three zeros 

W+{x) = -^^x{x^ - a^){x'^ + 6^). (44) 

In order to have the asymptotic behaviour of the function W^{x) in the 
vicinity of the points x = ±1 which corresponds to the case 2a, namely 
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— l/(a; — 1) and — l/(a; + 1), we choose the parameter a as follows 

This value of a gives a nonsingular behaviour of the potential energy in the 
vicinity of a; = ±1. In order to have the same derivatives in the points of 
zeros of the function Wj^{x\ namely W+(0) = W_^(±a), we put the following 
value for parameter h 

from which it follows that a^ > 3. Then choosing 

3(a^^ 1)2 
we obtain the superpotentials 



2e = iy;(0) = W',i±a) = ,,^, ^,, (47) 



a; / 2(a4 - 6a^ + 3) 2{a^ - 3)x^ _ 15{a^ - 3) \ 
^"^^ ~ 6 1^ (a2 - 1)2 + (a2 - 1)2 ~ (a2-3)x2 + 2a2J ^ ' 

X / 2(a^-6a2 + 3) 2(a2 - 3)a:2 15(a2 - 3) _ 12 \ 

'^^ 6^ (a2-l)2 ^ (a2-l)2 ^ (a2 - 3)x2 + 2a2 ^2-1^^^ 

and the nonsingular potential energy 

^ x^ / 2(a^-6a2 + 3) 2(a2 - 3)x' _ 15(a2 - 3) ^ 
"^ ^ 36 1^ (a2-l)2 ^ (a2-l)2 (a2 - 3)a:2 + 2a2 j ^ ' 

1 f2{a*-6a^ + 3) 6{a^ - 3)x\ 15{a^ - 3) 30ia^-3fx^ \ 



6\ (a2-l)2 (a2-l)2 ' (a2_ 3)3.2 + 2a2 ((a2 - 3)a:2 + 2a2)2y 

For this QES potential we know explicitly the wave functions of the ground 
and the third excited states 

/-/ N ^ // 2 o\ 2 , o 2x5/4 /-(a^ -3)a;^ + 2(a^-6a2 + 3)a;2\ 
V^o i^) = ^o((a' - 3)a;2 + 2a^f/^exp i-^ Uia^-lV ) '^ ^ 

xjx' - a2) /^ -(a2-3)x^ + 2(a^-6a2 + 3)x2 ^^ ^^^^ 
^^ ^") = ^- ((a2-3)x2 + 2a2)V4 ^^P ( T^(^T^Ty ) '^'^^ 

Note that at a2 = 3 this potential becomes exactly solvable and corresponds 
to the harmonic oscillator. 
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5 Conclusions 

We propose the SUSY method for constructing the QES potentials with 
arbitrary two known energy levels and corresponding wave functions. This is 
an extension of our SUSY method proposed in [|l^ |21| where QES potentials 
with the ground and first excited states were obtained. In the proposed 
method the function W^{x) plays the role of a generating function. Choosing 
different functions W^{x) we obtain different QES potentials V-{x). The 
two known wave functions V'(7(a^) and ip~{x) correspond to the eigenstates 
with zero energy and the energy e, respectively. The state numbers of these 
wave functions depend on the properties of function W+{x) as described 
in the cases 1, 2a, 2b which is summarized in the general case in section 
3. In section 4 we consider the explicit examples of QES potentials with 
the rational generating function W-^{x). These examples yield some new 
one-parametric QES potentials. At some special values of parameters these 
potentials become exactly solvable and correspond to the harmonic oscillator. 
We can consider various new generating functions W^{x) and obtain new 
QES potentials. One of the interesting possibilities is to consider the periodic 
functions W^{x) which lead to the periodic QES potentials. This problem 
will be the subject of a separate paper. 
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